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1. Introduction. 

The search for exact bound-state solutions of wave equations, relativistic 
or non- relativistic, has been an important research area in quantum me- 
chanics. However, over the past decades, problems involving the multidi- 
mensional Schrodinger equation have been addressed by many researchers. 
For instance, Bateman investigated the relationship between the hydrogen 
atom and a harmonic oscillator potential in arbitrary dimensions [1]. The 
JV-dimensional Kratzer-Fues potential was discussed by Oyewumi [2], while 
the modified Kratzer-Fues potential plus the ring shape potential in D- 
dimensions by the Nikiforov-Uvarov method has also been considered [3]. 
Very recently, some quantum mechanical properties of the pseudoharmonic 
oscillator were discussed by Agboola et al [4]. 

The Hulthen potential is one of the important short-range potentials 
which behaves like a Coulomb potential for small values of r and decreases 
exponentially for large values of r. The Hulthen potential has received ex- 
tensive study in both relativistic and non-relativistic quantum mechanics 
[5, 6, 7, 8]. Unfortunately, quantum mechanical equations with the Hulthen 
potential can be solved analytically only for the s-states [9, 10, 11]. How- 
ever, a number of methods have been used to find the bound-state energies 
numerically [10-16] and analytically [17, 18]. 

Recently, an extension of this study to a multidimensional Hulthen po- 
tential in a Klein-Gordon equation was presented by Saad [19]. The main 
idea of their investigation relies on using an approximation for the centrifu- 
gal term. 

Moreover, the simple conventional Nikiforov-Uvarov method [20], which 
received much interest, has been introduced for solving Schrodinger equation 
[6], Klein-Gordon [21], Dirac [22] equations. 

It is therefore the aim of this paper to present the approximate solutions 
of the Schrodinger equation with the Hulthen potential in D-dimensions for 
£ states using the conventional Nikiforov-Uvarov method. The paper 
is arranged as follows: the next section gives the Schrodinger equation in 
D-dimensions; section 3 gives a brief description of the Nikiforov-Uvarov 
method, while the next section presents the bound-state to the Hulthen 
potential. In section 5, some expectation values are calculated; and finally, 
the conclusions of the work are presented in section 6. 
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2.0 Schrodinger Equation in D-dimensional Spaces. 



Using the D-dimensional polar cordinate with polar variable r (hyperradius) 
and the angular momentum variable 6\, O2, #3, @d-2, 4> (hyper angle), the 
Laplacian operator in the polar coordinate r, 9±, 62, #3, Od-2, & of the R D 
is 



dr 



dr 



+ 



(1) 



where A.j-,(Q) is a partial differential operator on the unit sphere S ^ 1 
(Laplace-Betrami operator or the grand orbital operator) define analogously 
to a three-dimensional angular momentum[2] as A. 2 D (Q) = — £|> -(A?-) where 

^■ij = Xi ~5x~ ~ x i~§x ^ or a ^ Cartesian component Xi of the 
D-dimensional vector (x±, X2, xd)- 

The D-dimensional Schrodinger equation has the form [2,4]: 



h 2 

-^V 2 D * nem (r, ft) + V(r)V nem (r, Q) = E^ nim {r, ft) (2) 

where ji is the reduced mass and Ti is the Plank's constant. 

2.1 The Hyperradial Equation for the Hulthen Potential. 

The Hulthen potential is given as [7, 8, 9, 19, and 23]: 

— (XV 

where a is the screening parameter and Z is a constant which is identified 
with the atomic number when the potential is used for atomic phenomena. 
Inserting Eq.(3) into Eq.(2), we have 



dr 
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Seperating the variables Eq.(4) becomes 



,1-D 



d_ 

dr 



1(1 + D - 2) 



ft* 



2ju 



E + Za 



l-e~ 



E + Za- 



1 



^„£ m (r, 0) 

(4) 

Rndr) = 
(5) 
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and 

A 2 D (Wi m (n) + PYr(n) = o (6) 

where the seperation constant (5 is given as 

f3 = £(£ + N + 2); 1 = 0,1,2,... (7) 

Thus, we write Eq.(5) as follow 



E + Za- 



RntXr) = 0| 



(8) 

where R n e(r) is the hyper radial part of the wave function,E is the energy 
eigenvalue, I is the Orbital angular momentum quantum number. Eq.(8) is 
the D-dimensional hyper radial Schrodinger equation for the Hulthen po- 
tential. 

3. The Nikiforov-Uvarov Method. 

In this section, we give a brief description of the conventional Nikiforov- 
Uvarov method. A more detailed description of the method can be ob- 
tained the following references [20] .With an appropriate transformation s = 
s(r),the one dimensional Schrdinger equation can be reduced to a general- 
ized equation of hypergeometric type which can be written as follows: 

+ 44^00 + 4r~M s ) = ( 9 ) 

<J{S) a \ s ) 

Where a(s)and a(s) are polynomials, at most second-degree, and f(s)is 
at most a first-order polynomial. To find particular solution of Eq.(9) by 
separation of variables, if one deals with 



Eq.(9)becomes 
where 



1>(s) = Hs)yn(s), (10) 

a(s)y^ + T(s)y / n + Xy n = (11) 



w=*Wy (12) 

t{s) =f{s) + 2ir(s),T'(s) <0, (13) 
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, , a'(s)-f(s) I fa'(s) -f(s)\ 2 ~~ ~~ , 

<s)= y > 2 ^ ' ±\J[ - 2 — j -<7(a)+M«), (14) 

and 

\ = t + ir'(s). (15) 

The polynomial t(s) with the parameter s and prime factors show the dif- 
ferentials at first degree be negative. However,determination of parameter 
t is the essential point in the calculation of vr(s). It is simply defined by 
setting the discriminate of the square root to zero [20]. Therefore, one gets 
a general quadratic equation for t.The values of t can be used for calculation 
of energy eigenvalues using the following equation 

\ = t + tt'(s) = -nr'(s) - ^ ~ ^ a"(s). (16) 

Furthermore, the other part y n (s) of the wave function in Eq.(12) is the 
hypergeometric-type function whose polynomial solutions are given by Ro- 
drigues relation: 

- W)^ n(sMs)] (17) 

where B n is a normalizing constant and the weight function p(s) must satisfy 
the condition [20] 

{ap)' = rp. (18) 



4.0. Eigenvalues of the Hulthen potential in D-dimensions. 



We now want to obtain the solution to Eq.(8) using the Nikiforov-Uvarov 
method. If one defines 



(D-l) 

R n e(r) = r 2 U n e(r), 



(19) 



Eq.(8)becomes 
'2fi ( 



h 2 



E + Za 



1-e- 



(2£ + D - 1)(2£ + D - 3) 
4 r 2 



U ne (r) = 0. 
(20) 

Eq.(20) is similar to the one dimensional Schrodinger equation for the 
Hulthen potential, expect for the addition of the term ^(^+^~ 2 ) ;W hich is 
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well known as the centrifugal term. To solve Eq.(20),one can consider the 
approximation of the centrifugal term which is valid for a small value of a 



a 2 e~ ar 



(21) 



r 2 (l-e- Qr ) 2 ' 

With the use of Eq.(21) and the transformation s = e~ ar , Eq.(20) becomes 

U M+^^) U M+ [s{1 z s)]2 l(-e 2 S)s 2 + (2e 2 +5-~f)s-e 2 ]U ne (s) = 

(22) 

where 



2fiE 



2Zfi 



and 7 



{2i + D -1)(2£ + D -3) 



a 2 % 2 ' a% 21 ' 4 

By comparing Eqs.(9) and (22), we can define the following 



(23) 



f(s) = 1-s, a(s) = s(l-s) and a(s) = (-e 2 -5)s 2 + (2e 2 +6--/)s-e 2 

(24) 

Inserting these expressions into Eq.(14), we have 

7T(a) = -| ± \ + 4(e 2 + ^ - t)]s 2 - 4(2e 2 + 6 - 7 + t)« + 4e 2 (25) 

The constant parameter i can be found by the condition that the expression 
under the square root has a double zero, i.e. the discriminant is zero. Thus, 
the two possible functions for each t are given as 



' \ [(2e - VTT^f) s - 2e] for t = 7 - 5 + e^TT^y 



ir(r) = ± < 

v ' 2 



\ [(2e + yTT47>-2e] for t = 7 - 5 - e^/TT^i 

(26) 

However, for the polynomial r(s) = f (s) + 27r(s) to have a negative deriva- 
tive, we can select the physically valid solution to be 



7r(s) = -2"2 



(2e + v 7 ! + 47) s - 2e 



for i = 7 — 5 — e-y/1 + 47 such that 



r(s) = 1 - 2s - (2e + VTTI7) s - 2e 



(27) 



(28) 
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Also, by Eq.(16) we can define 

a = 7 - s - \{i + 2e) (1 + yrn^) 



n 



1 + 2e + n + Vl + 47 (29) 



for n = 1, 2, ... Thus, with the aid of Eqs.(23)and (29), after simple manipu- 
lations, we have the energy eigenvalue as 

2 



2*2 



En 



2jJL 



1 n(n + 2i + D - 2) + (7 - 5) 

2 + 2n + 21 + D - 1 



(30) 



where 7 and <5 are define in Eq.(23).We now give the following comment. 
Comments: Eq.(30) gives the energy spectrum of the Hulthen potential in 
D-dimensions. However, in order to check the validity of the approximation 
(21), we show that Eq.(30) can be reduced to a 3-dimensional case. That is 
when D = 3 and 7 = 1(1 + 1) and after a simple manipulation, we have 



E„, — — 



2fi 



n + e + i 



n + e + i 



-a 



(31) 



which is in agreement with previous works [6,24]. Also, we note that if we 
write Eq.(30) as follows 



E n — — — 



2\x 



a an(n + 21 + D - 2) + (07 
2 + 2n + 21 + D - 1 



2Zn\,-i 2 



and take the limit as a — > 0, we have 



En ~ 2h 2 



2Z 



1 2 



2n + 21 + D - 1 



(32) 



(33) 



which is the energy eigenvalue of the Coulombic potential [29]. Evidently, 
this follows from the fact that lim a ^o Vn(r) = — ^ . 



4.1. Eigenfunctions of the Hulthen potential in D-dimensions. 



In this section, we obtain the wave functions using the Nikiforov-Uvarov 
method. By substituting tt(s) and a(s) into Eq.(12), and solving the first 
order differential equation we have 

0(a) = s e (l -sf/ 2 , v = 2£ + D-l. (34) 

Also by Eq.(18), the weight function p(s) can be obtained as 

P (s) = s*(l-sy- 1 (35) 
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Substituting Eq.(35)into the Rodrigues relation (17), we have 



y n {s) = B n s 2e (l - s 



\1-V 



ds n 



Therefore, we can write the wave function U n i{s) as 

u ne (s) = c nS \\ - s y/ 2 pt' v - l) {i - 2S) 



(36) 



(37) 



where C n is the normalization constant, and we have used the definition of 
the Jacobi polynomials [26], given as 



(-If 



d" 



n!2"(l - s) a (l + s) b ds n 



(1 - s) a+n {l + s) b+n 



(38) 



To compute the normalization constant C n , it is easy to show with the use 
of Eq.(19) that 



/ \R ne (r)\ 2 r D - 1 dr= / 
Jo Jo 



U ni (r)\ 2 dr 



f 

Jo 



U nf (s)\ 2 — = 1 



as 



(39) 



where we have also used the substitution s = e ar . Putting Eq.(37) into 
Eq.(39) and using the following definition of the Jacobi polynomial[26] 



T(n + a+ 1) 
n!r(l + a) 



F\ I —n, a + b + n + 1; 1 + a; 



1 - s 



(40) 



we arrived at 



C 



r(2e + n+l) 



n!r(2e + 1) 



2 -i 



/ s 2e_1 (l - s) v 2 Fi (-n, 2e + v + n; 1 + 2e; s) ds = a 
Jo 

(41) 

where F is the hypergeometric function. Using the following series repre- 
sentation of the hypergeometric fucntion 



p F q (a 1 ,...,a p -c 1 ,...,c q ;s) = ^ — j-^-t (42) 



we have 



Cl 



r(2e + n + l) 



n!r(2e + 1) 



2 i-n)k(n + 2e + v) k (-n)j(w + 2e + v)j f 1 2e+k+j _ x _ _ I 

^ " (l + 2 e)fe fc! " (l + 2e) 3 jl Jo U s ^ s "-| 

(43) 
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Hence, by the definition of the Beta function,Eq.(43)becomes 



Cl 



r(2e + n + i; 



-, 2 n n 



n!r(2e + 1) 



E E y ~ nh ! n t?M ^ " B(2e+k+j, v+1) = a.\ 



k=0j=0 



(l + 2e) k k\ 



(l + 2e)jj\ 



(44) 

Using the relations B(x,y) = ^x+^y) an< ^ tne Pochhammer symbol 
(a) n = r r(o^ ' Eq.(44) can be written as 



r(2e + n + i; 



n 2 n 



n!r(2e + 1) 



\-n)k(2e)k(n + 2e + ^ (— n)j(2e + k)j(n + 2e + u)j _ a I 
^ (1 + 2e + + 2e) k k\ f^ Q (1 + 2e + v + A%(1 + 2e),j! " B(2e,v + l)| 



(45) 

Eq.(45) can be used to compute the normalization constants for n = 0, 1, 2, ... 
In particular for the gound state, n = 0, we have 



B(2e,u + 1) 



(46) 



5. Some Expectation Values for the Hulthen Potential in 
D-dimensons. 

We now calculate some expectation values of the Hulthen potential using 
the Hellmann-Feynmann theorem (HFT) [27,28]. Suppose the Hamiltonian H 
for a particular quantum system is a function of some parameters q, and let 
E n (q) and ^ n (q) be the eigenvalues and eigenfunctions of H{q) respectively, 
then the HFT states that 

°jM =mq)l ?m mq)) (47) 

The effective Hamiltonian of the hyperradial function is given as 

In order to calculate (r~ 2 ), we set q = I such that 

<*n(*)|^|*»(*)> = ^ + D ~ 2){r ' 2) (49) 



and 

dE n _ a 2 h 2 [(16<5 2 - (2n + 21 + £> - l) 4 ] 
~ (2n + It + D - l) 3 



(50) 
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Thus by HFT, we have 



(r- 2 ) = 



a 2 [(16<5 2 - (2n + 21 + D - 1) 



4 (21 + D -2){2n + 2i + D - l) 3 
Similarly, putting q = Z in Eq.(47),we obtain (V(r)) as 



<V(r)> 



2aZ 



(2n + 21 + D — 1) 



1 n(n + 21 + D - 2) + (7 - 5) 



+ 



2n + 21 + D - 1 



(51) 



(52) 



where 8 = 2Zfi/ah 2 and 7 = (2£+D 1 ^ 2e+D 3 ) , Finally, we note that the ex- 
pectation value (T),can also be obtained from the fact that E n = (T) + (V). 



6. Conclusions. 



In this paper, we obtained the solutions to the D-dimensional Schrodinger 
equation with the Hulthen potential within the framework of an approx- 
imation to the centrifugal term. In solving the hyperradial equation, the 
Nikiforov-Uvarov method was used, and the energy eigenvalues obtained was 
found to agree with the 3-dimensional case when D = 3. The eigenfuctions 
were also obtained in terms of the Jacobi polynomials and the normaliza- 
tion constants were also computed in terms of the hyper geometric functions. 
However, it is important to note that the approximation (21) is only valid 
for a small value of a; and as a — > 0, the results obtained approach those of 
the Coulombic potential. 

Moreover, some expectation values for the Hulthen potential in D-dimensions| 
were worked out using the Hellmann-Feynmann theorem. Finally, we note 
that the results obtained are parameter dependent. 
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